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ABSTRACT 

Several measurements of QSO-galaxy correlations have reported signals much larger 
than predictions of magnification by large-scale structure. We find that the expected 
signal depends stronly on the properties of the foreground galaxy population. On 
arcminute scales it can be either larger or smaller by a factor of two for different 
galaxy types in comparison with a linearly biased version of the mass distribution. 
Thus the resolution of some of the excess measurements may lie in examining the halo 
occupation properties of the galaxy population sampled by a given survey; this is also 
the primary information such measurements will provide. 

We use the halo model of clustering and simulations to predict the magnifica- 
tion induced cross-correlations and errors for forthcoming surveys. With the full Sloan 
Digital Sky Survey the statistical errors will be below 1 percent for the galaxy-galaxy 
correlations and significantly larger for QSO-galaxy correlations. Thus accurate con- 
straints on parameters of the galaxy halo occupation distribution can be obtained 
from small scale measurements and on the bias parameter from large scales. Since the 
lensing induced cross-correlation measures the first moment of the halo occupation 
number of galaxies, these measurements can provide the basis for interpreting galaxy 
clustering measurements which measure the second and higher order moments. 

Key words: cosmology: dark matter — cosmology: gravitational lensing — galaxies: 
clustering 



1 INTRODUCTION 

Gravitational lensing by large-scale structure along the line- 
of-sight can alter the observed number density of galaxies 
on the sky (Gunn 1967). Lensing increases the area of a 
given patch on the sky, thus diluting the number density. 
On the other hand, galaxies fainter than the limiting mag- 
nitude are brightened due to lensing and may be included 
in the sample, thus increasing the number density. The net 
effect, known as magnification bias, can go either way: it 
can lead to an enhancement or suppression of the observed 
number density of galaxies, depending on the slope a of the 
number counts of galaxies No{m) in a sample with limiting 
magnitude m, s = dlog No/dm. Magnification by amount 
H changes the number counts to (e.g. Broadhurst, Taylor & 
Peacock 1995), 

iV'(m) = iVo(m)^i^-^^-^ . (1) 

For the critical value s — 0.4, magnification does not change 
the number density; it leads to an excess for s > 0.4, and a 
deficit for s < 0.4. 



This magnification effect also leads to an excess corre- 
lation of QSOs and foreground galaxies associated with the 
lensing mass distribution (Canizares 1981; Schneider, Ehlers 
& Falco 1992; see also Keel 1982; Peacock 1982). It is mea- 
sured through the excess or deficit in the counts of galaxies 
around background QSOs. It is critical to ensure that the 
galaxies are not physically associated with the QSOs; this 
can be done by choosing appropriate cuts in the redshifts 
of the two populations. Several measurements of the QSO- 
galaxy cross-correlation have been made: significant excesses 
of galaxies around QSOs have been found, although the mag- 
nitude of the effect measured remains well in excess of the 
expectations of dark-matter theory with simple biasing, as 
discussed below (Bartelmann & Schneider 2001 give an ex- 
cellent review of the theory and observations). This has re- 
mained an unresolved problem, manifested in particular in 
measurements on angular scales of order an arcminute (see 
the discussion in Bemtez et al 2001). Calculations of the ex- 
pected amplitude of this effect have been done using both 
linear and non-linear dark-matter clustering (Bartelmann 
1995). 
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Several observational measurements of QSO-galaxy cor- 
relations have been made, beginning with the mid-1980's. 
While the results vary considerably, most studies with rsr 
dio loud QSOs have found correlations in excess of lens- 
ing predictions (Tyson 1986; Fugmann 1988,1990; Hammer 
& Le Fevre 1990; Hintzen et al. 1991; Drinkwater et al. 
1992; Thomas et al. 1995; Bartelmann & Schneider 1993b, 
1994; Bartsch, Schneider, & Bartelmann 1997; Bem'tez et 
al. 1995; Bem'tez & Martinez-Gonzalez 1995, 1997; Norman 
& Williams 1999; Norman & Impcy 1999, 2000). Although 
these results are qualitatively in agreement with the mag- 
nification bias effect, in most cases the amplitude of the 
correlation is much higher than that expected from grav- 
itational lensing models. Studies that have used optically 
selected QSOs have found both positive and null or nega- 
tive correlations. Webster et al (1988), Rodrigues- Williams 
& Hogan (1991), Williams & Irwin (1998) and Gaztanaga 
(2003) found positive correlations that are much stronger 
than the predictions from the magnification bias effect, 
while Boyle et al. (1988); Romani & Maoz (1992); Fcrreras, 
Benftez & Martinez-Gonzalez (1997) and Groom & Shanks 
(1999) found null or even negative correlations, which in 
some cases cannot be explained by the lensing hypothesis. 
Some of the measured cross-correlations are likely to be af- 
fected by the incompleteness of the QSO samples, and most 
have small samples of galaxies and QSOs with little infor- 
mation on the redshifts of the galaxy population. 

Even with these caveats, it is fair to say that there 
are severe discrepancies between measured QSO-galaxy cor- 
relations and theoretical predictions. Two kinds of resolu- 
tions have been discussed in the literature: observational 
effects, such as bias in the selection of QSO samples, phys- 
ical associations of QSOs and galaxies, and dust obscu- 
ration; and improved theoretical modeling of the lensing 
(see the discussions in Benftez, Sanz & Martinez-Gonzalez 
2001 and Bartelmann &: Schneider 2001). On the theory 
side, various nonlinear effects have been modeled. Bartel- 
mann & Schneider (1993) and Bartelmann (1995) provided 
the framework for theoretical modeling in terms of large- 
scale structure parameters. Dolag & Bartelmann (1997) and 
Sanz et al (1997) made magnification predictions including 
nonlinear gravitational evolution, as done by Jain & Seljak 
(1997) for the shear, while WilHams (2001) and Menard et 
al (2002) studied corrections to the weak lensing approxi- 
mations. However, severe discrepancies remain; in a careful 
analysis, Benftez, Sanz & Martinez- Gonzalez (2001) find a 
signal that is a factor of a few larger than model predictions 
on scales of 1 — 5'. 

In this paper we model the galaxy distribution using the 
halo model of clustering (e.g. Cooray & Sheth 2002). We 
compute QSO-galaxy and galaxy-galaxy cross-correlation 
using this model for the galaxies and compare it to the linear 
bieis approach used so far in the literature. We seek to ex- 
plain some of the discrepancies between measurements and 
lensing predictions with our galaxy model. Under the lens- 
ing hypothesis, the QSO-galaxy cross-correlation is due to 
the cross-correlation of galaxy number density with magni- 
fication. This in turn can be expressed as the projection of 
the galaxy-mass cross-correlation. We only modify the next, 
final step, of expressing the gala^xy-mass cross-correlation as 
a linear bias factor times the mass-mass correlation. Instead 
we use the halo model of clustering to populate halos of given 



mass with galaxies, and then compute the cross-correlation 
consistently. One might expect that on small enough scales, 
^ 1 Mpc, the results would differ significantly from a linear 
bias model because the occupation of halos with galaxies 
is complicated and depends strongly on galaxy type. Our 
approach to galaxy-mass correlations is similar to that of 
Seljak (2000), but there are differences in the galaxy model 
used. 

We aim our predictions at forthcoming survey data, pri- 
marily from the Sloan Digital Sky Survey (SDSS), and the 
CFHT Legacy survey. We make predictions for the QSO- 
galaxy and the galaxy-galaxy cross-correlation. Section 1 
contains the formalism for the cross-correlation calculation. 
Section 2 details the halo model and the prescription for 
assigning galaxies to halos. The results are presented in Sec- 
tion 4 and we conclude in Section 5. 



2 FORMALISM 

Angular correlations can be expressed as projections of the 
3-dimensional power spectrum using Limber's equation. We 
follow the convention and derivations of Moessner & Jain 
(1998). Let rti {(/)) be the number density of foreground galax- 
ies with mean redshift (zi), observed in the direction (f> in 
the sky, and n2 (0) that of the sample with a higher mean 
redshift {z^} > The angular cross-correlation function 
is then defined as: 

uj{e) = {5ni{4>)5n2{4>')) , (2) 

where Sni{4>) = {ni{4>) — ni)/ni, with fii the average number 
density of the ith sample. The measured Srii is the sum of 
fiuctuations due to the true clustering of galaxies Snf, and 
due to magnification bias Sn^. 

The fluctuations on the sky due to intrinsic clustering 
are a projection of the galaxy fiuctuations along the line- 
of-sight, weighted with the radial distribution W{x) of the 
galaxies, 

SnU4>) = n dxWi{x)S''{r{x)4>,a) , (3) 
Jo 

where a denotes the expansion scale factor. The comoving 
radial coordinate x s-^d the comoving angular diameter dis- 
tance r(x) are explicitly defined in e.g. Jain & Seljak (1997). 
Xh is the distance to the horizon. 

In the weak lensing limit the magnification is /i = 1 + 
2k, where the convergence k is a weighted projection of the 
density field along the line-of-sight (see equation 5 below; see 
also Menard et al 2002 for tests of the weak lensing limit). 
Equation 1 then gives the magnification induced fluctuations 
in number counts as 

5<(0) = 5(s-O.4)Ki(0), (4) 
with the convergence k given by 

= n dx 5i(x)^^^ , (5) 

where the radial weight function g{x) can be expressed in 
terms of r{x)\ for ^ delta- function distribution of source 
galaxies at xs it is g{x) = r{x)r{xs - x)/r{xs)- 

The cross-correlation oj{ff) is then given by the sum of 
four terms. In the case where the background and foreground 
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samples have no overlap, the cross-correlation is dominated 
by (Moessner & Jain 1998), 



3n^(2.5s2 - l^Tv 



I 



X / dkkPgm{x,k)Jo[kr{x)9] , 



(6) 



where Pgm(x>^) is the galaxy-mass cross-power spectrum. 

Note that the expression for the mean tangential shear 
around a foreground galaxy is identical to the above equa- 
tion except for the replacement of Jo by J2 (this follows 
from the relation 7t(^) = —1/2 dK{9)/dln9). The difference 
in Bessel functions that appear for the shear and magni- 
fication can be useful in combining information from the 
two measurements. They weight length scales in different 
ways because the magnification measures probe the local 
projected density while the shear measures probe the inte- 
grated mass inside an aperture. Thus with comparable signal 
to noise measurements from the two probes, one can gain not 
only an imrovement in overall sensitivity, but also probe a 
larger range of length scales. 



3 MODELS FOR GALAXY-DARK MATTER 
CORRELATIONS 

We express galaxy-mass cross-correlations by developing a 
halo model based description of the cross-power spectrum, 
Pgm(Xi k), introduced above. Pgm can be modeled in at least 
three ways: (i) Take it to be a bias factor times the matter- 
matter power spectrum; (ii) Assume there is a one-to-one 
correspondence between galaxies and halos, and that the 
galaxy sits at the centre of its parent halo; (iii) Use the halo 
model (Seljak 2000; Peacock & Smith 2000; Scoccimarro et 
al. 2001; see the recent review by Cooray & Sheth 2002), 
which allows for the possibility that the number of galaxies 
in a halo may vary from halo to halo, and allows for the pos- 
sibility that the distribution of galaxies around the parent 
halo centre may depend on galaxy type. 

In the halo model, all dark matter particles and galaxies 
inhabit halos. This means that 



Pgm(fc) — 



Pgm{k) -\- Pgm(^) 



Pgm(fe) 



-I 



dmn(m) ^ (-^gail^) |wgg^i(fe|m)| 



Pin(fc) j 



P '^gal 

dm n{m) — b{m) u{k\m) 



X / dmn(m) -^^^^i^i^^ 6(m) Ugai(/c|m) 



dmn(m) m 



ngal 



dm n{m) {Nga,\\m) . 



(7) 



Here n(m) is the comoving number density of halos which 
have mass m, (A*'gai|'m) is the average number of galaxies in 
halos of mass m, h(m) is the bias parameter, and u{k\m) 
and ttgai(fc|m) are the Fourier transforms of the density run 
of dark matter and galaxies around the centres of the halos 
they inhabit. The models assume that the density run can 



be well parametrized by allowing for a depedence on halo 
mass: 



u{k\m) 



1 J dr Anr'^ p{r\m) sm{kr)/kr 



(27r)3 



J dr 47rr^ p(r|m) 



(8) 



where the factor of (27r)^'^ defines our Fourier transform 
convention, and the integral in the denominator defines 
the mass m of the halo. The denominator normalizes 
(27r)^w(A;|m) so that it is unity at small k; for all profiles 
of interest, it decreases to zero as k increases, although the 
decrease need not be monotoruc. 

Option (ii) — one and only galaxy per halo, and the 
galaxy sits at the halo centre — corresponds to setting 
A'gai(m) = 1 and (27r)'^'Ugai(fc|m) = 1 for all m (we use 
Ngaiim) to deuote (Afgai|m) for simplicity). This is a rea- 
sonable model for the Large Red Galaxy (LRG) sample pro- 
vided by the SDSS; results for this model are shown in Fig- 
ure 4. 

Option (iii) requires additional assumptions to fully 
specify the galaxy distribution. For example, one might ex- 
pect the distribution of galaxies around the halo centre to 
depend on galaxy type. In what follows, we will assume for 
simplicity that the galaxies in a halo trace the dark matter: 
therefore, we set Mgai(fc|m) = u{k\m). Sheth et al. (2001) 
show that this is a reasonable approximation. (See Scran- 
ton 2001 for more discussion of how Wgai might depend on 
galaxy type). 

It is probably accurate to assume that one of the A^gai 
galaxies in a halo sits at the halo centre - this is almost 
certainly a good assumption if there is only one galgixy in a 
halo. To account for this, one assumes that the other galax- 
ies get the usual density profile denoted by u, whereas the 
central galaxy is simply unity. To see the effect of this on the 
two-halo term, we must average both pieces of the two halo 
term over p(n|m), the probability of having n galaxies in a 
halo of mass m, with the requirement that n > 0. With this 
averaging, we will replace A'gai(m)ugai(fe|m) by an effective 
number of galaxies, iVefi, which is a function of k for given 
m. It is given by (Cooray & Sheth 2002) 

Nes{k\m) = ^ [1 + {n - 1) u{k\m)] p{n\m) 

+ ^ (n - l|m) u{k\m) + u{k\m) p(0|m). (9) 

7l>0 

The first term can be written as [1 — p{0\rn)], allowing us to 
express the above equation as 

iVeff(A;|m) = [l-p(0|m)] [l-Wgai(/c|m)]-h(n|m) Mgai(/c|m).(10) 

Since both factors in the first term are positive, there is 
an enhancement in power from always placing one galaxy 
at the halo centre. Since u{k\m) decreases eis k increases, 

the enhancement in power is largest on small scales (large 
k). In sufficiently massive halos one might expect to have 
many galaxies, and so p(0|m) <g; 1. In this limit, the ex- 
pression above becomes 1 — u{k\m) + (n\m) u{k\ni) = 1 -|- 
(n — l|m) u{k\m). On the other hand, if most halos have no 
galaxies, then p(l|m) is probably much larger than all other 
p{n\m) with n> 2. Then the leading order term in the sum 
above is p{l\m). Since (n|m) = ^np(n|m) w p(l|m), we 
have that Nes{k\m) « (n|m). In this limit, only a fraction 
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(n|m) <C 1 of the halos contain a galaxy, and the galaxy sits 
at the halo centre, so there is no factor of u. 

The contribution of the galaxy counts to the one halo 
term of the galaxy-mass correlation function is similar. Us- 
ing the expressions above yields 



-fgm(fe) 



-Pgm(fc) 



-I 



dmn{m) — \u{k\m) 



Wgal 



Plin(fc) 



dmn(m) — b{m) u{k\m) 



dmn{m) b{m) 



N^s{k\m) 



(11) 



If the run of galaxies around the halo centre is not the 
same as of the dark matter, then one simply uses Mgai in- 
stead of u in A^eff. If the two-halo term usually does not 
dominate the power on small scales (this is almost always 
a good approximation), it is reasonable to ignore the en- 
hancement in power associated with the central galaxy, and 
to simply sot Ncti{k\m) ~ {n\m) Ugi^i{k\Tn) ~ {n\m). The 
one-halo term requires knowledge of p(0|m). Since p(0|m) 
is usually unknown, some authors (e.g., Seljak 2000) inter- 
polate between the two limits discussed earlier by setting 
A'eff = (n|m) u{k\m) if {n\m) > 1, and A''eff — (njm) if 
(n|m) < 1. It is worth noting that the one halo term of 
the galaxy-galaxy correlation function is only slightly more 
complicated, as shown in Cooray & Sheth (2002). 



3.1 Halo model details 

The halo model is spocifiod by the mass function n(m), bias 
6(m) and the halo profile u{k\m). For the first two functions 
we use (e.g. Cooray & Sheth 2002) 



n(m, z) dm 

h{m) 
where v 



m V 



A{p) (1 - 
qv- 



(qu) 



')(£)""-(-f) 



1 + 



1 ^ 2p/6,,{z) 



4o(2) 1-^(91^)^'' 



Sljz) 



(12) 



andp w 0.3, A(p) = [1-^2-^ r(l/2-p)/V7r]"^ « 0.3222, and 
q w 0.75 (Sheth & Tormen 1999). Here 4o(«) is the critical 

density required for spherical collapse at z, extrapolated to 
the present time using linear theory, and 



<7 (m) = 



An 
(2^ 



Jo 



dk 

^A;^PLi„(fc) H/'(fcJ?o), 



(13) 



where W{x) = (3/x^)[sin(a;) 



sm a; — xcosia; 



and Ro 



(3m/47rp)^^^. That is to say, a{m) is the rms value of the 
initial fluctuation field when it is smoothed with a tophat 
filter of comoving size Rq, extrapolated using linear theory 
to the present time. If p = 1/2 and q — 1, then n(m) is 
the same as that first written down by Press & Schechter 
(1974), and 6(m) is the same as that given by Cole & Kaiser 
(1989) and Mo & White (1996). 

In addition, we will assume that the halo profiles have 
the form given by Navarro, Frenk & White (1996), trun- 
cated at the virial radius rvir which is defined by requiring 
that m = 47r/3rvirp Avir. For spatially flat universes with 



fio = (1, 0.3) and A = 1 - Avir = (178, 340). The Fourier 
transform of the density run around a halo of mass m is 



u{k\m) = f{c) 



+ cos 



i/t(^Si[/t(l-|-c)] -Si(K)j 
/t(^Ci[/t(l-|-c)]-Ci(/t) 



i(kc) 



(14) 



where /(c) — l/[ln(l + c) — c/(l -I- c)], k = krvir/c, 
Si(x) = dtsm{t)/t is the sine integral and Ci(x) = 
— dt cos(t)/t is the cosine integral function. The con- 
centration parameter of the halos depends on halo mass; we 
use the parametrization of this dependence given by Bullock 
et al. (2001): 



c(m) w 9 



(-)■ 



(15) 



3.2 Galaxy Model Details 

Our model for the gala^xy distribution is taken from 
the GIF N-body simulations of the ACDM model, 
coupled to a semi-analytic galaxy formation model 
(Kauffmann et al. 1999). Catalogs of galaxy positions, 
apjjarcnt and absolute magnitudes, colors and star- 
formation rates, at a range of redshifts, are available at 
http://www.mpa-garching.mpg.de/GIF/. From the cata- 
logs at redshifts z = 0.06, 0.13, 0.27, 0.35, 0.42, and 0.52 
we selected apparent magnitude limited samples which sat- 
isfied 17 < r* < 21. At each redshift wo then computed the 
first and second factorial moments of the number of galaxies 
as a function of halo mass. The parameters for the first mo- 
ment are given in Table 1 . These galaxies also form a surface 
defined by the first and second factorial moments as a func- 
tion of halo mass and galaxy absolute magnitude which we 
can also integrate over given our apparent magnitude limits 
and a redshift (Scranton 2002). 

In the GIF models, the mean number of galaxies in 
a halo (n|m) evolves little out to z ~ 0.5 if the galaxies 
have the same fixed rest-frame luminosity. Thus at fixed 
luminosity the evolution of galaxy clustering is driven by the 
evolution of the halo population. However, this means that 
at fixed apparent magnitude, the case in which wo arc most 
interested, {n\m) evolves more strongly. This is because only 
the most luminous galaxies from the higher redshift samples 
satisfy the apparent magnitude Umit, and, in halos which 
host only one galaxy, there is a tight correlation between 
halo mass and galaxy luminosity. In principle, fits to (n|r7i) 
at fixed luminosity, along with the assumption that these 
relations do not evolve, allow one to estimate how {n\m) 
depends on redshift for the apparent magnitude limits we 
apply. In practice, we chose the more tedious but accurate 
method of making redshift dependent fits to the apparent 
magnitude (n\m) relation in the GIF model. The results are 
shown in Figure 1. 



3.3 Lens and source redshift distributions 

Magnification bias with galaxy-galaxy lousing should be de- 
tected with SDSS data and forthcoming datasets such as 
the CFHT Legacy survey. The SDSS photometric survey 
will provide of order 1 million foreground galaxies and 100 
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Table 1. Parameters which specify the mean number of galax- 
ies in halos as a function of halo mass: Ngi^i{m) s (m/Mo)" + 
j4exp 1^— j4o (/ogio(m) — Mg)^j . This Table gives the parameter 
values for the qso-galaxy case, with the foreground galaxies se- 
lected in the magnitude range 17 < r' < 22. 



Table 4. Parameters for redshift distribution functions of galaxies 
and QSOs of the form dN/dz = exp[-(2/2o)'']- 



2 


log Mo/h~^MQ 


a 


A 


Ao 


Mb /log (/i 


0.06 


12.30 


1.05 


0.98 


11.56 


11.55 


0.13 


12.47 


1.03 


0.96 


0.62 


11.42 


0.27 


12.74 


0.94 


0.87 


1.65 


12.10 


0.35 


12.82 


0.79 


0.84 


7.41 


12.25 


0.42 


13.11 


0.76 


0.91 


7.3 


12.45 


0.52 


13.67 


0.84 


1.03 


10.07 


12.68 



2 1 


log Mo/h~^MQ 


a 


A 


^0 


Ms/log {h-^Mo) 


0.06 


12.26 


0.9 


1.0 


2.0 


11.75 


0.13 


12.48 


0.9 


1.0 


4.0 


11.75 


0.27 


13.0 


0.9 


1.0 


4.0 


12.75 


0.35 


13.48 


0.9 


1.0 


6.0 


12.75 


0.42 


14.0 


0.9 


1.0 


8.0 


12.75 


0.52 


11.:-! 


D.I) 


l).2 


8.0 


12.75 



million background galaxies with photometric rcdshifts. The 
expected signal is sensitive to the mean redshifts of the fore- 
ground and source populations, but only at the 10% level or 
less to the shape of the distribution. Given that there is 
some freedom in the selection of the two populations with 
photometric redshifts, we have chosen simple distributions 
for our model predictions. For a CFHT Legacy survey kind 
of dataset, we choose the foreground population to be uni- 
formly distributed over 0.2 < z < 0.5 and the sources at 
z = 1. For predictions for the SDSS survey, we choose a 
much closer foreground population with 0.1 < z < 0.2, and 
the sources aX, z = 0.4. We have checked that for a realis- 
tic, broad source redshift distribution with the same mean 
redshift, the signal would be lower by about 5%. 

Our expected distributions of galaxies and QSOs from 
an SDSS-like dataset for measurement of the QSO magnifi- 
cation bias are much more constrained than for the galaxy- 
galaxy lensirig case. First, we expect that the efficiency of 
photometric selection of QSOs, while high compared to pre- 
vious work, will force the use of spectroscopically-confirmed 



Table 3. iVgai(m) as in Table 2, but for "red" foreground galaxies, 
defined by the criterion g' — r' > 0.65. The functional form is 
iVgai(m) s (m/Mo)"exp(-Mi/m). 



z 


log Mo/h^'^MQ 


a 


log Mi/h-'^MQ 


0.06 


11.70 


1.2 


10.0 


0.15 


11.70 


1.0 


11.7 


0.26 


12.18 


0.8 


12.0 


0.35 


12.18 


0.5 


12.0 


0.42 


12.85 


0.5 


12.0 


0.52 


12.85 


0.5 


12.0 



Object 


A 


a 


20 


b 


Galaxies 


1 


1.56 


0.296 


1.76 


Low-2 QSOs 


253.28 


2.58 


1.67 


13.44 


High-2 QSOs 


6.5 X 10-* 


9.13 


3.35 


18.37 



Table 2. A'gai(m) as in Tabic 1, but for the galaxy-galaxy case, 
with the foreground galaxies in the apparent magnitude range 
17 < r' < 21. 



QSOs. This will result in approximately 10® QSOs, down 

three orders of magnitude from the number of galaxies in 
the photometric sample. Hence we will need to maximize the 
number of objects in our distributions rather than restrict- 
ing them to narrower, easier-to-model bands in redshift. For 
galaxies, this means going to a broad cut in apparent mag- 
nitude: 17 < r' < 22, for example. Using the luminosity 
function from the CN0C2 survey (Lin et al 1999), we can 
translate such an apparent magnitude cut in the SDSS fil- 
ters into a redshift distribution (Dodelson et al 2002) well 
fit by a simple functional form: 



cx z" exp f— (z/zo) 



(16) 



(IN 

dz 

The QSO sample distribution is more complicated. Due 
to a combination of the QSO luminosity function and the 
evolution of QSO colors in the SDSS filters with redshift, the 
redshift distriubtion of QSOs from the main SDSS sample 
requires two functions of the form given in Equation 16, with 
an appropriate amplitude (^4) for each to properly charac- 
terize the relative abundance of QSOs below z = 2 and those 
above. Table 4 gives the parameters for the galaxy and QSO 
redshift distributions, with the normalization set to A = 1 
for the galaxy distribution. In addition to these smootii dis- 
tributions, we impose a strict redshift cut oi z < 0.8 for the 
galaxy distribution and 2; > 1 for the QSO distribution to 
avoid intrinsic clustering between the two samples. 



4 RESULTS 

Figures 2-4 show the predicted cross- correlation for differ- 
ent foreground- background populations. We use equation 6, 
with an f2m = 0.3, Q,\ = 0.7, F = 0.21 cosmological model 
with as = 0.8 and the galaxy-mass power spectrum deter- 
mined using the halo model ingredients described in Section 
3. While erg = 0.8 is close to the most recent determinations 
from the CMB and other datasets, there is uncertainty at 
about the 10% level in this parameter (e.g. Wang et al 2002; 
Spergel et al 2003; Verde et al 2003). Our results, especially 
for the higher redshift lens galaxies, are sensitive to ag: for 
higher erg, the lens galaxies are less biased and vice versa. 
Note that we have left out the (2.5s — 1) factor which de- 
pends on galaxy type. It will modulate both the amplitude 
and possibly the sign of our plotted cross-correlation. 

We consider galaxy-galaxy correlations in Figure 2 
for redshift distributions appropriate to the CFHT Legacy 
survey in the upper panels and the SDSS survey in the 
lower panel. The right panels show the ratio of the cross- 
correlation to the Peaeock-Dodds prediction. On large scales 
this reduces to the usual bias parameter, but on small scales 
it must be interpreted more carefully in terms of proper- 
ties of the halo occupation distribution. The results show 
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that the cross-correlations arc strongly dependent on galaxy- 
type and angular scale. This arises from the modeling of the 
"red" and "blue" galaxy sub-samples that make up the full 
foreground sample. The galaxy samples are defined by the 
Nga,i{m) models given in Tables 2 and 3 for the full sam- 
ple and for red galEixies (the blue galaxies are given by the 
difference of the two). Similar trends are seen in Figure 3, 
which uses the galaxy distribution in Table 4 and a mean 
redhisft z = 1.75 for the quasars. 

It is interesting that for red galaxies, lowering the red- 
shifts of the source and lens population does not change the 
signal significantly (comparing the upper and lower panels 
of Figure 2). This is a consequence of how little Afgai(m) 
evolves between z — 0.5 and z — for these galaxies. The 
net change in lensing signal as one lowers the redshifts is a 
trade off between two opposing effects: a reduction due to 
smaller lensing path length, and an increase because a given 
angle corresponds to a smaller length scale at the lens red- 
shift. With a simplified power law model for the galaxy-mass 
power spectrum one can see that the net signal can increase 
with decreasing redshift. Let Pgm(fe, a) oc k~°'a'^, and let 
the lenses and sources be at a single redshift such that the 
distance to sources is twice the distance to the lenses. Then 
for an Einstein-de Sitter model equation 6 gives the analytic 
scaling 

w{e) oc a-t^^" e'-". (17) 

Thus for a = 1 (corresponding to a two-point correlation 
function with logarithmic slope of —2), one can see that 

/3 > 1 would load to an increasing signal with decreas- 
ing lens redshift. The choice a = 1 and /3 = 1 is in good 
agreement with the behavior of the red galaxy sample, but 
not of the blue galaxies for which the cross-correlation is 
lower for lower redshifts of the lenses and sources. The dif- 
ferences between red and blue galaxies arise primarily from 
the difference in the Ngai{m) relation for them (Sheth et 
al 2001; Scranton 2002). On larger scales, corresponding to 
r > IMpc, the difference between predictions for red and 
blue galaxies decreases. In the linear regime they are both 
expected to follow the growth rate of mass fluctuations, so 
the large-scale differences are primarily due to differences in 
a. 

We estimate the contribution of Poisson errors and sam- 
ple variance on the measured u}(0). Scranton (2002) discusses 
the relative contributions of these errors and of the Gaus- 
sian and non-Gaussian terms in the covariance. Over the 
scales we have considered, the Poisson contribution domi- 
nates for a large survey like the SDSS. For parameters of 
the full SDSS survey, roughly 10 million galaxies in each 
sample and an area of 10,000 square degrees, the statisti- 
cal errors are tiny. The errors are of order 1% of the signal 
and have been multiplied by a factor of 10 in Figure 2 to be 
visible at all. Clearly the signal to noise is high enough for 
such a dataset that estimates for sub-samples of galaxies by 
type and varying redshift bin are possible, thus allowing for 
the possibility of measuring their halo occupation proper- 
ties in detail. The signal to noise can be scaled to other sur- 
vey parameters relatively easily on small scales, where shot 
noise dominates so that the scaling is oj/Slj oc ^nin2%/0, 
where ni and 712 are the number densities of the two sam- 
ples, and n is the survey area. Figure 3 shows that the errors 



expected from quasar-galaxy correlation measurements are 
significantly higher, due to the smaller number of quasars. 
Thus if systematic errors can be controlled, magnification 
bicis is more efFectly measured using galaxy-galaxy correla- 
tions. 

The two major source of systematic errors are photo- 
metric calibration and errors in the photometric redshifts 
that are not well characterized or highly correlated with 
determination of the galaxy type. The latter leads to con- 
tamination of the samples: e.g. a foreground galaxy may 
be assigned to the background sample in the galaxy-galaxy 
case, or may in reality be physically associated with the 
quasar. In either case the auto-correlation signal, which is 
much stronger than the lensing signal, will get mixed in. 
Moessner & Jain (1998) estimate the required accuracy for 
photometric redshifts for a desired accuracy in the cross- 
correlation. Given the miniscule statistical errors, charac- 
terizing the systematic errors accurately will be paramount 
for any interpretation of results from either the CFHT or 
SDSS surveys. 

For the quasar-galeixy u){d) shown in Figure 3 the error 
bars are larger because of the smaller number of quasars 
compared to background galaxies. Even so, they arc at the 
1% level and would allow for a definitive measurement from 
the SDSS provided systematic errors are under control. 

Figure 4 shows the cross-correlation for halos of differ- 
ent masses. This is not an observable, but it is helpful to 
sec the wide variation in the signal as a function of halo 
mass since galaxies are modelled as occupying these halos. 
Clearly small changes in the iVgai(m) relation will change 
the expected signal since it will populate halos of different 
masses with galaxies differently. The dot-dashed curve shows 
the prediction for halos of mass lO"Af0, which h avc a very 
large amplitude at small scales. This is relevant for the SDSS 
since it will produce a sample of Large Red Galaxies, which 
arc expected to occupy the centers of massive halos. A more 
quantitative connection is beyond the scope of this paper. 

The redshift distribution of galaxies we have used has 
not been chosen to be consistent with the models for 
■^gai("^). This amounts to requiring from survey measure- 
ments that some (typically small) fraction of galaxies be 
discarded before estimating the cross-correlation to ensure 
consistency. 



5 DISCUSSION 

We have used the halo model of clustering to estimate 

the cross-correlation induced by magnification bias between 
samples of galaxies and quasars at different redshifts. Our 
focus has been on the effect of the model for galaxies used 
on the predicted signal. We used fits to the GIF simula- 
tions for iVgai(m), the mean number of galaxies for halos 
of mass m, to find the number of galaxies as a function 
of redshift, magnitude and galaxy type. We find that the 
predicted cross-correlation is very sensitive to these param- 
eters. Its amplitude varies on arcminute scales is a factor 
of 2 larger or smaller compared to an unbiased population, 
depending on the redshift range and galaxy type of the fore- 
ground population. The reason for this wide range is that 
on small scales the relation of the galaxy distribution to 
the mass is complex and cannot be described by a constant 
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bias factor of order unity. Further, with increasing redshift 
a given apparent magnitude limit corresponds to brighter 
absolute magnitudes and therefore the selection of galaxies 
with higher m/m*, which tend to be more strongly clus- 
tered. Our results on arcminutc scales therefore differ from 
previous studies which used linear bias models. 

The state of measurements of qso-galaxy correlations 
are summarized in Bem'tez et al (2001) and Guimaraes et 
al (2001). Clearly some of the measurements are affected by 
incomplete sampling and other observational systematics. 
Others find excess correlations on large scales, which are 
still difficult to reconcile with models. But most anomalous 
measurements arc on arcminutc scales whore we find great 
sensitivity to the redshift distribution and galaxy type used 
for the measurement. Thus variations of a factor of 2 to even 
10 in extreme cases can be explained. In existing data, the 
galaxy samples used are difficult to characterize given the 
lack of rcdhift or color information. The data discussed in 
Benftez & Martinez-Gonzalez (1995) show that the sample 
with red galaxies does have higher signal - this is consistent 
with our results, but it is difficult to be quantitative given 
the partial information available for the data. 

In view of the massive improvement in the available 
sample expected from the SDSS and other forthcoming sur- 
veys such as the CFHT Legacy survey, we have chosen to 
simply provide model predictions for these surveys rather 
than attempt a detailed analysis of past measurements. Our 
models can bo extended significantly, to include more de- 
tailed galEixy types and galaxy models that are checked for 
consistency with SDSS measurements of the galaxy auto- 
correlations. 

The statistical errors expected from the SDSS survey 
are very small, about 1% of the signal for the full sample 
from galaxy-galaxy correlations. We have focused on scales 
below 20 arcminutes, where the signal to noise is high. On 
these scales the data from the SDSS can be split by galaxy 
type and redshift bins to make a detailed study of galaxy 
clustering in relation to the dark matter distribution. On 
large-scales the interpretation is simple, giving the bias pa- 
rameter of desired galaxy sub-samples as a function of red- 
shift. 

What can we learn from sub-Mpc scale measurements of 
the galaxy-mass correlation? In the halo model, the 1-halo 
term is a linear measure of (the product of) A'gai(m) and 
other parameters: the halo mass function, halo mass pro- 
file and galaxy profile. Scranton (2002) analyzed the accu- 
racy with which different parameters describing the galaxy 
distribution can be constrained using angular correlations 
from the photometric SDSS survey. Our results show that 
valuable additional information can be gained using mea- 
surements of magnification bias. In particular, the lensing 
induced cross-correlation is a measure of projected galaxy- 
mass correlations: it is thus linear in the parameters of the 
galaxy distribution such as the mean halo occupation num- 
ber of galaxies A'gai(m). Galaxy clustering measurements 
probe the second moment of this distribution, and therefore 
one is required to make an assumption of how the second 
moment scales with the mean (on large scales galaxy cluster- 
ing measurements do constrain the first moment through the 
two-halo term, but only averaged over a broad mass range). 
The magnification bias measurements considered here would 
constrain iVgai (m) directly and thus provide the basis for in- 



terpreting the higher order moments from galaxy clustering 
measurements. Further, the lensing measurement would help 
break the degeneracies between the various parameters. 

Scranton (2002) gives the accuracy with which more 
than 10 parameter combinations will be constrained by 
SDSS measurements of the angular clustering of galEixies 
binned by photometric redshift. For the basic parameters 
of iVgai(m) in Table 1, he typically finds accuracies of well 
below 1%. For the lensing induced cross-correlation, the sig- 
nal is significantly smaller and for the dominant parameters 
in A/gai (m) such as a one can expect better than 10% level 
accuracy from the lensing measurements alone. The main 
power however will be in combining the lensing and galaxy 
clustering measurements. We leave this exercise for future 
work. 

Our formalism is similar to that of SeljaJj (2000) and 
Guzik & Seljak (2001) for galaxy-galaxy lensing. These au- 
thors have used measurements from the SDSS to constrain 
the mass profiles of large galaxy lialos and the contribu- 
tion of group halos (Guzik & Seljak 2002; see also McKay 
et al 2001). While in principal magnification bias measures 
a closely related cross-correlation (replacing the tangential 
shear with the convergence), in practice the measurements 
have been relevant on different scales. Galaxy-galaxy lensing 
has been measured so far on relatively small scales, the best 
signal being on 100 — 200kpc. It remains to be seen from the 
completed SDSS survey and CFHT surveys how the relevant 
systematics play out for magnification and shear measures. 
It will be of great interest to analyze the measurements 
jointly to constrain galaxy-mass correlations. Similarly the 
large scale bias parameter inferred from shear surveys (Van 
Waerbeke 1998; Hoekstra et al 2002) can be directly com- 
pared with measurements of magnification bias. This will 
be an extremely useful cross-check since the errors in both 
measurements are likely to be dominated by systematic er- 
rors. 

We have used the weak lousing approximation in this 
paper, setting /j, = 1 + 2k. There are corrections to this re- 
lation if K or 7 are of order unity (e.g. Menaxd et al 2002). 
Correlations using the fully nonlinear magnification rela- 
tion fi — 1/[(1 — k)^ — 7'^] can be computed using the halo 
model approach developed in Takada & Jain (2003a, b). It 
leads to higher amplitudes for magnification bias on subar- 
cminute scales, especially for source redshifts ^ 1 (Takada 
& Hamana 2003); it will therefore be important to incor- 
porate the nonlinear effects for deep surveys. Further, we 
have treated galaxies as points within halos, rather than 
subclumps with mass profiles. Using the formalism of Sheth 
& Jain (2002) we have estimated the effect and found it to 
be small for the range of scales studied in this paper. 

We thank N. Benftez, A. Connolly, R. Scoccimarro, A. 
Szalay and M. Takada for helpful discussions. We thank the 
referee, Brice Menard, whose comments helped improve and 
correct the paper. Some of the work presented here was be- 
gun during a summer 2002 workshop at the Aspen Center 
for Physics. BJ is supported by NASA grants NAGS- 10923, 
NAG5-10924 and a Keck foundation grant. 
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Figure 1. The mean number of galaxies as a function of tialo mass, Ng^i(m), is shown for a range of redshifts. The symbols show the 
measurements from GIF simulations, and the solid curves shows our fits in Table 1. 
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Figure 2. (i) The left panels show the angular cross correlation computed using the linear mass power spectrum (lower dotted curves), 
the nonlinear mass power spectrum (upper dotted), and three galaxy models (blue, total, red from bottom to top). The foreground 
galaxies are taken to lie between 0.2 < 2 < 0.5 (0.1 < 2 < 0.2) and the sources at 2 = 1 (0.4) for the upper (lower) panels. The error 
bars have been multiplied by 10 (see discussion in the text), (ii) The right panels show the ratio of the galaxy cross-correlations to the 
PD prediction. On large scales it is reasonable to interpret this ratio as a bias parameter, but on small scales it is important to consider 
the details of the halo occupation distribution of galaxies. 
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Figure 3. The angular cross correlation as in Figure 2, but for the qso-galaxy parameters appropriate for the SDSS survey. 
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Figure 4. The angular cross correlation predictions for halos of masses 10^*^, 10^^, 10^^, IO^^Mq (bottom to top), compared to the linear 
and PD predictions (solid curves). 
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